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^ ; ABSTRACT 

The cross sections for the production of B s and B* mesons through the frag- 
' mentation process in pp collisions were calculated. The input parameters were 

determined from the measured B s and B* production rates in e + e~ annihilation. 

The transverse momentum distributions da/dp? are given for yfs = 1.8 TeV, the 
\ Fermilab Tevatron energy. 
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At LEP the DELPHI Collaboration has measured the probability for a b antiquark 
| to hadronize into a weakly decaying strange B meson.El Assuming the hadronization is 

dominated by the fragmentation of the heavy b antiquark, we determined the fragmen- 
tation functions Di^ Bs (z, jj) and Di__> B * (z,pS) from the measured total probability for 
a b antiquark to hadronize into the lowest strange-beauty states B s (}So) and £>*( 3 Si).i 
Here z is the fraction of the b antiquark momentum carried by the B s or B* at the 
scale fi. The momentum distributions of the B s and B* mesons produced in e + e~ 
annihilation at the energy of the Z mass were then predicted. 
| We then extended our calculations to the production of the B s and B* states in 

^ pp annihilation at the Tevatronja using the fragmentation functions previously deter- 

minedJ 

In the parton model the cross section for the production by fragmentation of 
strange B mesons in proton-antiproton annihilation involves three main factors: the 
structure functions of the initial partons in the proton and antiproton, the subprocess 
in which the initial partons produce a particular parton in the final state, and the 
fragmentation of this final parton into the strange B meson. The transverse momentum 
distribution of the B q meson is then of the form 



I dxidxj dz fi/pix^fj/pixj)-^- (i(xi)j(xj) k(— )X, fi) D k _ Bs {z; //) . 
i,j,k J apr \ z J 

(1) 



*Talk presented at the International Conference on Fronties of Physics: Looking to the 21 st Century, 
Shantou, P.R. China, 5-9 August 1995. 
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Here fi/ p {xi) and fj/p{xj) are the structure functions of the initial partons % and j carry- 
ing fractions of the total momenta Xi and Xj in the proton and antiproton, respectively. 
The production of parton k with momentum p T /z is described by the hard subprocess 
cross section a at the scale fx. And Dk~,B s { z 'i /-0 i s the fragmentation function for the 
parton k to yield the meson B s with momentum fraction z at the scale fx. The factor- 
ization scale fx is of the order of the transverse momentum of the fragmenting parton. 
The soft physics is contained in the structure functions, while the parton subprocess 
is a hard process and can be reliably calculated in perturbation theory. We have also 
treated the fragmentation of b into B s as if it was also a hard process and perturbation 
theory was applicable even though the s quark is not very heavy. The main dependence 
on the s quark mass is hi the normalization of the fragmentation functions, which was 
empirically determined.H 

Clearly, there is a formula very similar to Eq. ([!]) for the production of the B*, 
which we need not write explicitly. Because the radiative decay B* —>■ B s + 7 is so fast 
and the B* — B s mass difference is so small, the transverse momentum of B* is hardly 
affected in the decay B* —>■ B s + 7, and thus contributes to the inclusive production 
of the B s . We will, therefore, present each individually in the figures showing da /dp? 
and <j(pt > £>™ in ) and also include their sum in a table, for convenience. 

In the fragmentation process k — > B s we have included gluons as well as b anti- 
quarks. Although the direct g — > B s fragmentation process does not occur until order 
there is a significant contribution coming from the evolution of the b antiquark frag- 
mentation functions from the heavy quark mass to the collider energy scale Q, which is 
of order log(Q/mf,) due to the splitting g — > b. This has been shown to be significant 
for the production of B c mesonso and it was also included. 

For the parton distribution functions we used the most recent CTEQ version 
3,0 in which we chose the leading-order fit, since our calculation is also a leading- 
order one. For the subprocesses we used the tree-level cross sections, to be consistent, 
since the parton fragmentation functions were calculated only to leading order. For the 
production of b antiquarks we included the processes gg — > bb, gb — > gb, and qq — > bb, 
while for the production of gluons g we included the processes gg — > gg, gq(q) — > gq(q), 
and qq — * gg, where q denotes any of the quarks, u, d, s, c, b. 

For the running strong coupling constant a s we used the simple one-loop result 
evolved from its value at \i = m^: 

a 'W = - , 33-2^ — ; — x ; 7 a \ ■ ( 2 ) 

Here rif is the number of active flavors at the scale \x and we chose a s {mz) = 0.1 17.il 
The scale \x superficially appears to have been invented only to separate the 
production of B s mesons into structure functions, subprocess cross sections, and frag- 
mentation functions. In principle, this is so and the production is independent of the 
choice of \i. But in practice, the results do depend on \i, since only if all the factors are 
calculated to all orders in a s will the dependence on /i cancel. However, for example, 
the fragmentation functions were calculated only to leading order and, consequently, 
the results will depend on fx. We investigated the sensitivity of the results to the choice 
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of [i by varying fj, from /xr/2 to 2/xr, where [Lr = yp^(parton) + ml is our central 
choice of the scale \i. 

To obtain the fragmentation functions at the scale /z we numerically integrated 
the Altarelli-Parisi evolution equations^ from the scale /io, which is of the order of the 
6-quark mass. As boundary conditions in these calculations we used the fragmentation 
functions Di_> Bs (z, /io) and Di_> B * (z, /io) at the scale /i calculated using perturba- 
tive QCD.i The induced gluon fragmentation functions D g ^ Ba>B *(z, /i) are of order 
al log(/i/ra&) and become important at large values of the scale \i relative to the b anti- 
quark fragmentation functions Di^b s ,b*{ z i A*)' which are of order a^(fi). The boundary 
conditions for the gluon fragmentation functions are D 9 ^b s (z, /•*) = D g ^ B *(z, A-0 = 
for jj, < 2(mb + m s ), the threshold for producing B s or B* mesons from a gluon. 

We chose the b quark mass to be m h = 5 GeV and the strange-quark mass 
parameter m s was determineda using the same initial fragmentation functions from 
the experimental valued of the probability ff for a b antiquark to fragment into weakly 
decaying strange-beauty mesons: 

f: = f 1 dz [D- b _ Bs (z, /i ) + D- b _ Bt (z, /i )l . (3) 

Since the total probability for the b antiquark to fragment into a B meson is indepen- 
dent of scale, the initial scale /io in Eq. (|3]), which is of the order of the 6-quark mass,, 
was chosen to be /io = rrib+2m s . Then using the measured value ff = 0.19±0.06±0.08S 
and a s (m z ) = 0.117a we then obtain m s = 298 + 47- 23 MeV. 

The evolution equations were numerically integrated and the fragmentation func- 



tions at the scale /i# = ypy(parton) + mf were combined with the structure functions 
and parton cross sections for the subprocesses to obtain the cross section, Eq. (0). In 
Fig. |I]we show the cross section da/dpx for both B s and B*. Specifically, we calculated 
the transverse momentum distributions of B s and B* mesons produced in pp collisions 
at the Tevatron energy y/s = 1.8 TeV. We assumed a cut-off on the transverse momen- 
tum of px > 5 GeV/c 2 and considered only the rapidity range \y\ < 1. To investigate the 
sensitivity of these results to the scale \i we have also included the results for /i = 2/xr 
and fi = Hr/2. When the scale \i is less than fi = nib + 2m s , which only happens for 
the case of \i = Hr/2, we chose the larger of (//, A*o)- From Fig. |l| it is clear that the 
choice of scale \x is not critical for the transverse momentum distribution; in fact, the 
variation over the range /xr/2 < ji < 2//r is comparable to the current discrepancies 
between the measured and calculated b production cross sections. As one might expect, 
the ( 3 «Si) B* cross section is larger than the ( 1 5 , ) B s cross section at all pt, however, 
their ratio is not precisely the naive prediction 3:1, but is about 20% smaller. 

In Fig. we show the total cross section for the production of B s and B* mesons 
with transverse momentum above a minimum value p™ m . As in Fig. |l| only the range 
Pt > 5 GeV and \y\ < 1 was considered. The sensitivity to the scale \i was investigated, 
as before, by considering /i = Hr/2 and \i = 2/i R , and the choice of this scale is clearly 
not critical. Figure |] provides useful estimates of the production rates of B s and B* 
mesons at the Tevatron due to the fragmentation process. 
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p T (B a ) (GeV) 



Fig. 1. The transverse momentum distributions of B s and B* mesons for pr > 5 GeV, 
\y\ < 1, and fj, = /j.r/2, fiR, and 2//r at ^fs = 1.8 TeV. 

For convenience, we have also summarized these cross section estimates, including 
the sum a(B s ) + cr(£>*), in Table 1. 

Table 1. Total B s and B* cross sections in nb versus p™ m for px > 5 GeV and \y\ < 1 at 
-y/s = 1.8 TeV. Variation with the scale /i are shown. 
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The production rates and transverse momentum distributions for B s and B* me- 
son production in pp collisions presented here assume that fragmentation is the domi- 
nant process and, therefore, are probably underestimates. At large enough transverse 
momentum the fragmentation process should dominate, as it falls off more slowly, even 
though it is only a part of the full order a A s contribution. A comparative study of the 
relative importance of the various contributions in the full order a A s calculation is now 
being carried outp which will clarify the range of validity of the fragmentation approx- 
imation. And the comparison with forthcoming data from the Fermilab Tevatron will 
be very instructive. 



4 




These calculations were carried out in collaboration with Kingman Cheung and 
are reported in greater detailed in references [2] and [3]. This work was supported 
by the U. S. Department of Energy, Division of High Energy Physics, under Grant 
DE-FG02-91-ER40684 and DE-FG03-93ER40757. 
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